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PREFACE 


This paper is a study of some of the relationships between geomag- 
netic activity, zonal index, and meridional index. The superposed- 
epoch method is used to try to find relationships between mean zonal 
index and various types of geomagnetic activities. The same method 
is used to try to find relationships between mean meridional index and 
the various types of geomagnetic activities. 

This work was conducted at the United States Naval Postgraduate 
School, Monterey, California during the winter am spring of 1954 in 
partial fulfillment of the requiremerts for the degree of Master of 
Science in Aerology. 

The author wishes to thank Professor F. L. Martin, Department of 
Aerology, United States Naval Postgraduate School for his guidance. A 
varticular debt of gratitude is due to Professors A. B. Mewborn and 
C. L. Perry for their assistance in connection with the statistical 


analysis. 
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TABLE OF SYMBOLS & ABBREVIATIONS 


D= declination (or "variation") of the magnetic field 
F = total intensity of the earth's field 

H = horizontal component of the earth's field 

B= vertical component of the earth's field 

I = inclination of the earth's field 


~ North = South component of the earth's field 


East - West component of the earth's field 


<~"« 
it 


107? oersted 


K = geomagnetic index expressing local features such as the sys- 
tematic diurnal variation in geomagnetic activity 

K = "standardized" indices freed, insofar as possible from local 
features 

Kp= average of the K, indices from eleven (at present) magnetic 
observations 

MI= meridional index from 80° W to 0 


Ange Slat. change in height (in feet) from 80° W to the first trough or 
Sow R 
ridge at the particular latitude at which MI is computed 


[DhTn of = change in height (in feet) from the last trough or ridge 
‘ (before reaching 0 longitude) to zero longitude 


Ah§ = change in height (in feet) from preceding trough or ridge 
to first downstream ridge or t rough 

Do= day of increase in the twenty-four hour Kp index by nine or 
more units and a drop of this index within seventy-two hours 
by nine or more units 


Subscript i =, th gay of observation where Oth day is the "key" day, 
18% day is "key" day 4 1, etc. 


th 


Zi-= zonal index of an observation in the sample on the i— day 

ni- total number of observations in the sample on the it? day 

Zzmi= = 724i - mean value of the zonal index for the re day 
ni 
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S212 - zg + = anim of the mean of the sample on the 
nic yth day 
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- unbiased estimate of the variance of the popy lation based 


on the variances of the samples of all the i— days 
Zmi - Zmo x standard normal viriate using 9S mo as an ap- 
O mo proximation to the standard deviation of the 
population 
Zmi - Zmo - standard normal variate usingOm as unbiased 
m estimate of standard deviation of the popula-. 
tion 


+ 
meridional index of an observation in the sample on the i” da 


FZ MIi_ = mean value of the meridional index for the it? day 
et 
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CHAPTER I 


INTRODUCTION 


In the past six years many investigations of the effects of geomag- 


“netic activity on various meteorological elements have been carried out. 


- 


Most of the investigations have been concerned with the study of geomag- 
netic activity and its effect on pressure. Several of these have lead 
to interesting conclusions. 

Duell and Duel1° have studied the variation of surface pressure in 
northwest Europe during the few days following Seite icAlty disturbed 
days and also following geomagnetically quiet days. Their investigation 
led to the result that surface pressure variations during these periods 
were linked with events in the upper atmosphere. Craigt, at the Harvard 
College Observatory, felt the great need in the field was for an objec- 
tive statistical analysis with an aim toward providing informtion as to 
the probability that the Duell's results did demonstrate some physical 
relationship. His statistical analysis led him to state "It therefore 
seems reasonable to reject the hypothesis that geomagnetic activity and 
surface pressure are unrelated," 

Kaciak and Langwe11® have studied the relation between geomagnetic 
activity and pressure variations aloft. They did not restrict their in- 
vestigation to winter months and years of low relative sunspot numbers 
as did the Duells' and Craig. Their results indicated no relationship 
between geomagnetically disturbed days and pressure aloft. There did 
seem to be a relationship between geomagnetically quiet days and pres- 


sure aloft - especially at and north of 60° N. 


















At this point an investigation as to the relations between peomag- 
netic activity and circulation indices was both natural and d’ sired. io 
restrictions were made as to sunspot numbers. In some of the investi- 
gations no restrictions as to seasons were made. In others, only the 
winter season was investigited. The "winter" season was chosen as the 


months January, February, M.rch, April, November, and December. 








CAAPTER. II 


GEOMAGNETIC ACTIVITY 


To describe the earth's magnetic field completely at a point, the 
direction and intensity of the field at that point must be known. Figure I 
shows schematically the elements which define the marpnetic field at a given 
point. The three magnetic elements D (declination or "variation"), H(hori- 
zontal component) and Z(vertical component) are normally recorded at ob- 
servatories. The total intensity, F, of the earth's field can be deter- 
mined by computation from the relationship 

FaHSeé I (261). 
The inclination I can be measured by a dip circle or earth inductor. It 
is also easily calculated from the relation 
Tan I = —s (28>) 
The north component X and the east component Y are computed by the equations 


X=HCosD (2,3) 


Y=-H SinD (2.4) 


The unit of magnetic intensity is the oersted. In most places on 
the surface of the earth the intensity of the earth's magnetic field is 
less‘than one oersted. It is convenient in the taking of measurements 
to use the unit gamma (J) which is 0.00001 oersted. 

The total intensity F of the earth's magnetic field is not constant 
in time. The intensity varies and according to the rates of these var- 
dations are classified as secular or transient variations. The secular 
variations produce perceptible effects over periods of decades or cen- 


turies. Their origin is unknown. The transient variations have periods 
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measured in days, hours, or even minutes. There are three types 
of transient variations namely: (a) solar (b) lunar (c) disturbance 
variations. We are particularly interested in the disturbance varia- 
tions am their relation to the circulation indices. 

. The intensity of these disturbances are measured every three hours 
by a so called K-index. The range for K® 9 is 500 gamma (one famma 2 
eO0001 oersted). The numerical character figure of the K index is 
determined from the magnétogram trace every three hours. The K index 
expresses the degree of irregulzrities in the systematic diurnal var- 
lation of geomagnetism which includes certain local features. There 
was a need for an abstract of all the individual K-indices in order 
to express world wide features in the time-chanpes of feomagnetic 
disturbance. In 1949 this need was met by the establishment of a 
new index Kp designed to measure planetary variations in geomarnetic 
activity. This index was based on "standardized" K, indices freed, 
insofar as possible, from local features. «At present,eleven obser- 
vatories have conversion tables assigning a Kg index to their res- 
pective K indices. K is the average’ of the e¥even K, indices. The 
eleven observatories are: Lerwick, Meanook, Sitka, Hskaldemuir, 

Rude Skov, Agincourt, Vingst, witteveen, Abinger, Cheltenham, and 
Amberly. 

The indices K and Ky have been defined as measures of the dis- 
turbance exhibited by magnetograms. Actually these indices are de- 
signed to measure the varying intensity of solar particle radiations 
P by their geomagnetic effects. These solar particle radiations show 
their effects mostly in the auroral zones (polar regions). The eleven 


selected observatories lie within forty-two degrees of the respective 























poles of the earth's magnetic axis. According to J. Bartels, of the 
University of Géttingen, eruptive processes in the so called M regions 
of the sun generate these particle invasions. These particle invasions 
are not well correlated with sunspots, bright spots, or solar flares. 
Ionospheric storms (characterized by a considerable decrease of the 
electron density in the F layer) are well correlated with geomapnetic 


disturbances, 





CuArTsER Lil 


MSTHOD CF COLLECTION OF DATA 


All zonal profile valucs were obtained throuch project Arowa from 
the National Weather Records Center at Asheville, North Carolina. The 
records contained the mean geostrophic westerly wind for 5° latitude 
belts from 10% to 70N for the four octants and the hemisphere. The 
octant descriptions were given as; 

i. LOOK to LEO 

2. 17O0l to 9OW 

3. 8OW to 0 

4h. OE to 90E 
Values are in meters por second. Octant three (80°W to 0°) was selec- 
ted as the zone of investigation due to the frequency of blocks in this 
area and also due to the fact that some of the investigation involved 
the use of data from the observatory at Cheltenham, Maryland. This oo- 
servatory is on the western edge of octant three. 

Meridijonal indices were calculated from the Weather Bureau Daily 
Series 500 mb synoptic weather maps. The indices were calculated for 
octant three by starting at £0°Wand proceeding to 0°. As the flow 
changed from south to north or north to south the height cnange, in 
feet, from trouch to ridge or ridge to trough was recorded. These 
height changes, regardless of sign, as recorded from 80°W to 0° were 


a. 


added ‘tie 
HT = [Obes ORY + [OAS] t]OKE tee +[OhRa ec] (3.1) 


+* This sum was the fiscure actually used in tne computations of means 
and standard deviation. 





To change these values to meters per second multiply the sum of the 


heights by 


si. - 3048 
(“r-)“80(Gos BY 60 x 1853. ~—S—«(3-2) 


where g = 9.8 meters per second per second; f = 2#Sin @ = 1.454 X (Gn 
Sin %.; and # = latitude. Latitudes selected for both zonal and meri- 
dional indices were 50°N and 65°N. 

All geomagnetic data was taken from the Journal of Geophysical 
Research, The sum of the Kp index (value of the intensity of distur- 
bances for a twenty-four hour period) for each of the five most quiet 
days per month was recorded for the entire year 1949. This data was 
also broken down so that only the winter months were considered. The 
days so recorded were designated as "quiet key" days. The twenty- 
four hour Kp index for each of the five most disturbed days per month 
was recorded for the years 1949 and 1950. These years were also 
broken down so that only the winter months were considered. The diys 
so recorded were designated as "disturbed key" days. 

In selecting "quiet" or "disturbed key" days in this manner it 
was found that many of the "key" days fell on successive days. Since 
there was a possibility that this could introduce overlapping influ- 
ences, other methods of selecting "key" days were used. The first of 
these methods was the choice of days of s.c. (sudden commencement ) 
or me (small initial impulse followed by main impulse) as "key" days. 
These days were selected from Cheltenham, Maryland's records of prin- 
cipal magnetic storms for the years 1949, 1950, and 1951. One s.c. 
occurred on December 27, 1951 at Cheltenham that was not included in 


the sample. This was due to the fact that at the time the sample 











was being taken, data for this month and year were not available. A 


second method of selecting well separated "key" days was the selection 






of the day of increase in the twenty-four hour Kp index by nine or more 
units and a drop of this index within seventy-two hours by nine or more 
units. In general this will give a "key" day commencing on the day of 
increase in particle invasions after they had been relatively low. 
These days were selected for the winters of 1949, 1950, and 1951 and on 
the Plates are labeled "Dp". A third and fourth method of selecting 
"key" days were used. They involved using the twenty-four hour Kp in- 
dex for the single most disturbed day per month and the twenty-four hour 
Kp index for the single most quiet day per month. These days were chose 


for each month for the years 1949, 1950, and 1951. 





CHAPTER IV 


APPLICATION OF SUPERPOS“D EPOCH METHOD 


In order to see the relation, if one exists, between geomagnetic 
activity and either of the circulation indices all other causes of 
variation in the circulation indices must be eliminated or held constant. 
This means that the statisticil approach is the most natural and best 
suited for the investigation. Such a method has been widely used in the 
past and will be used here. It is known as the Chree "superposed-epoch" 
Method. 

It is necessary first to select certain well defined days of geo- 
magnetic activity called "key" days and secom try to determine if the 
geomagnetic activity of these days does or does not affect the particu- 
lar parameters we are investigating. Zonal and meridional indices were 
chosen as the parameters to be iny@@taested. Tne values of these para- 
meters over a particular octant (80°W to 0°) and at particular latitudes 
(50°N and 65°N) are recorded for each "key" day as well as several days 
following the "key" day. Py averaging the zonal and meridional indices 
over all the key days, a mean value for each type of index is obtained 
for that day. Likewise a mean value of the indices is obtained for each 
of the days following the key day. In this way we get a typical behavior 
pattern of the indices for the periods during and after the key day. 
This pattern should represent the effect of geomagnetic activity on the 
zonal indices provided all other causes of index variations are elimi- 
nated. 


There are three major types of index variations: diurnal, annual, 


10 





and interdiurnal. The interdiurnal variations are due to variations in 
the synoptic situation and are assumed to be eliminated due to the size 

of the sample. The diurnal index changes have been eliminated by using 

a single value computed for the same period each day. The annual varia- 
tion, in some cases, has been partially eliminated or minimized by using 
the same amount of data each month. In other cases only the winter months 
are considered and therefore an effect characteristic of the winter scason 
is incorporated in the results. 

If the subscript i designates the ith day of observation where i #0 
designates the "key-day", i#1 designates "key-day" plus one; i=2 desig- 
nates "key-day" plus two; etc. and Zi = zonal index of an observation 
in the sample on the ia day then the mean value of the zonal index for 
the pal day is 

mi = 2h (4,1) 


where Ni = total number of observations in the sample on the pas day. 


sett 
The variance of the mean of the semple on the gth day is estimated by 


ad cd b 


Let us éssume that, although the expected value of the population 
of means for each of the various days (i = 0, 1, 2, 3, 4, 6, 8) may not 
be the same for each value of i, the variance of the population of means 
does not depend on i. Under this assumption we can obtain an unbiased 
estimate of the variance of the population of means based on the variances 


of the samrles of all the es days. 


tt 
Mood (8) page 133 


dB 





It can be written as 


Sin =No ia od N — Gs 6m +N ie i Sn a (4.3) 
mr Re Rat ee, . a 
where Ni enuals the number of observations in the sample of the im day. 
The means and variances of the meridional index can be written in the 
same way using MI for meridional index. 

Let us assume that we have an gen day sample of individual zonal 
index (or meridional imiex) values of sufficient size that the sample 
means may be considered to be normally distributed. If the null hypo- 
thesis that the mean of the zonal index (or meridional index) for the 
key day is equal to the mean of the index for any of the days following 
the key day is made, this essentially states that there is no relation- 
ship (obtainable with this method and data) between geomagnetic activity 
and zonal index (or meridional index), If however, testing shows that 
Zmi = Zmo =1.966m, then there are only five chances out of a hundred 
that a deviation of the sample means as large or larger than this could 
occur due to chance alone. If Zmi - Zmo=1.966m it is due to the 
probability, with 95% confidence, that the samples are not randomly 
selected and the null hypothesis must be rejected. In view of the 
‘manner of selection of key days this implies that the deviation may be 
attributed to geomagnetic activity. 

Using the standard normal variate tr = ae and the table dn 
Hoel*, page 243, we find the probabilities of a deviation less than 


certain values of the standard normal variate is as shown in Table l. 


3 
Hoel (4) Page 130 
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TABLE 1 


The application of the Chree Method is shown in Tables 2 through 
ll, which present values of the standard normal variate following the 


particular choice of key day noted. In all tables an asterisk indicates 


a value of the standard normal variate in excess of 1.966m. 
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CHAPTER V 


DISCUSSION OF RESULTS AND CONCLUSIONS 


The results of the Tables 2 through 11 are depicted in graphical 
form in Plates I through IV. The standard Ak timesoH O mo was used 
first as a measure of the standard deviation of the mean of the popu- 
lation. Black dashed lines were drawn at distances £O mo, 220 mo so 
that it would be a matter of easy inspection to determine the days 
which give deviations in the neigiborhood of £20mo. A better esti- 
mate of the standard deviation of the population of means,Om of 
ecuation 1.3, was then devised, and the appropriatet On, + 20m boun- 
daries delineated on the plates in red dashed lines. Tests of signif- 
icance relative to the 20m limits are probably to be preferred although 
in some cases, there is little difference betweenOm andO mo. 

Plates I and II show the variations of zonal indices following vari- 
ous selections of "key" days and Plates III and IV show the variations 
of meridional indices following the same selections of "key" days as 
on Plates I and II. On Plate I, sections (3), (5), and (7) (essentially 
cases of disturbed geomagnetic activity) show variations at 50° N on the 
fourth day in excess of 20m. On Plate II, section (2) (most disturbed 
days) also shows a fourth day variation at 50°N in excess of 20m, In 
addition, Plate II, section (3) (dealing with another well marked type 
of disturbance) shows a similar trend at 50 N but the amount of varia- 
tion on the fourth pou is not quite at the 20m level. These results 
indicate that the value of the mean zonal index at 50N should be con- 


siderably greater on the fourth day following geomagnetic disturbances 


20 





than on the "key" day. Note the general increase in the value of the 
mean zonal index occurring between "disturbed key" day and "key" day 
plus four. This increasing trend, occurring as it does in all tests 
at SO°N involving disturbed days appears to indicate a significant 
development in the general circulation. 

Plate I, sections (1) and (2), showing variations of zonal indices 
at SO°N following "quiet key" days in 1949 has no variation greater than 
20m except on the eighth day. This implies that the value of the mean 
zonal index will be greater on the eighth day following quiet days than 
on the "key" day. This trend toward an increase on the eighth day can 
be seen independently from Plate II section (1), although not at the 20m 
level. The increase between roughly the third and eighth days following 
"quiet key" day is similar in slope to that occurring just after "dis- 
turbed key" day. This seems to indicate that the zonal index variation 
eight days after "quiet key" day could be due to the occurrence of a 
geomagnetically disturbed day at about quiet day plus three since on the 
long term average these should be 3 days apart in a 30 day month. 

For tests at 65°N Plate I section (8) and Plate II section (4) imply 
that the value of the mean zonal index drops to a minimum on the fourth 
day following geomagnetic disturbances. It is interesting to note that 
the minimum at 65°N and the maximum at 50°N occur on the same day fol- 
lowing geomagnetic disturbances, namely the fourth day. Meteorologi- 
cally this seems to imply that the circulation is altered in the direc- 
tion of a single strong-peaked zonal profile, corresponding to the 


classical "high index" stage (cf. for example Chapter II of "Forecasting 
in Middle Latitudes" by Riehl’). 


2) 





The meridional indices on Plates III and IV do not seem to have as 
repular or well pronounced variations as the zonal indices. Tests giving 
variations in excess of 20m or a variation very close to this value are 
shown in Plate III section (4) and Plate IV sections (1) and (3). Sec- 
tions (1) and (3), Plate IV, do suggest a maximum decrease in the mean 
meridional index on the third or fourth day following the "key" day. 
This appears to be an independent substantiation of the results concern- 
ing zonal index since synoptic experience indicates that the classical 
"high index" is accompanied by weak meridional index. 

This study indicates that a relationship exists between disturbed 
geomagnetic activity and the peneral circulation. The relationship 
indicated is that following a geomagnetically disturbed day one might, 
on the average, expect a trend toward the classical "high index" pattern 
at the 500 millibar level from 0° longitude to 80°W (i.e. octant three) 
and along 50°N latitude. | 

For further work along this line, fruitful results might also be 
derived from the consideration of (1) a different selection of "key" 
days (2) an investigation of the behavior of circulation indices on 
days before "key" days (3) a study of geomagnetic activity as related 


to CIT Weather Types or changes in Weather Types. 
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